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In a previous work [Phys. Rev. A 85, 022502 (2012)] we calculated, with the use of the Gaussian expansion 
method for few-body systems, the energy levels and spatial structures of the 4 He trimer and tetramer ground 
and excited states using the LM2M2 potential, which has a very strong short-range repulsion. In this work, 
we calculate the same quantities using the current most accurate 4 He- 4 He potential [M. Przybytek et al, Phys. 
Rev. Lett. 104, 183003 (2010)] that includes the adiabatic, relativistic, QED and residual retardation correc- 
tions. Contributions of the corrections to the tetramer ground-(excited-)state energy —573.90 (—132.70) mK 
are respectively -4.13 (-1.52) mK, +9.37 (+3.48) mK, -1.20 (-0.46) mK and +0.16 (+0.07) mK. Further 
including other realistic 4 He potentials, we calculated the binding energies of the trimer and tetramer ground 
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and excited states, B% ,B\ ,B\ and B\ , respectively. We found that the four kinds of the binding energies 
for the different potentials exhibit perfect linear correlations between any two of them over the range of binding 
energies relevant for 4 He atoms (namely, six types of the generalized Tjon lines are observed). The dimerlike- 
pair model for 4 He clusters, proposed in the previous work, predicts a simple interaction-independent relation 

„(i) „(0) 

+ |, which precisely explains the correlation between the tetramer excited-state energy and the 



B 2 



trimer ground-state energy, with B-2 being the dimer binding energy. 



I. INTRODUCTION 

The bosonic J = + three and four 4 He atom systems, 
which are very weakly bound under the 4 He- 4 He potential 
with an extremely strong repulsive core followed by the van 
der Waals attraction, are known to be suitable for studying the 
Efimov effect and the universality in the systems interacting 
with large scattering length |Q]-|5|]. 

In a previous paper [6], referred to as paper I in the follow- 
ing, we presented state-of-the-art four-body calculations for 
the 4 He tetramer ground- and excited-state binding energies 
and structural properties using a realistic 4 He potential called 
LM2M2 1 7], which has a very strong short-range repulsion. 
At the same time, our three-body calculation reproduced all 
the well known results for the 4 He trimer. We took the Gaus- 
sian expansion method (GEM) for ab initio variational calcu- 
lations of few-body systems nSJ- TTOfl . The total wave function 
is expanded in terms of totally symmetrized few-body Gaus- 
sian basis functions, ranging from very compact to very dif- 
fuse with the Gaussian ranges in geometric sequences. 

The method is suitable for describing the short-range cor- 
relations (without a priori assumption of any two-body corre- 
lation function) and the long-range asymptotic behavior (see 
the review papers [ 10-13] for many applications of the GEM) . 
As a result, we found in paper I that precisely the same shapes 
of the short-range correlation (r,j < 4A) in the dimer appear 
in the ground and excited states of the trimer and tetramer 
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and that the wave functions of the very weakly-bound ex- 
cited states of the trimer and the tetramer reproduce the correct 
asymptotic behavior up to up to ~ 1000A. 

Recently, Przybytek et al. UM proposed a 4 He pair poten- 
tial that is currently most accurate. Such an accurate 4 He 
potential is of importance, according to Ref. [14], in sev- 
eral branches of science, for example, in metrology (ther- 
modynamics standards) Il5l4l7ll . helium-nanodroplet spec- 
troscopy J[18l LLfaL and low-temperature condensed matter 
physics J20I1 as well as in the study of the unusually large 
and very weakly bound states of the 4 He clusters. The 
potential of Ref. [ 14] includes, in addition to the standard 
Born-Oppenheimer (BO) potential, various post-BO contri- 
butions. The main contributions are (i) the adiabatic correc- 
tions resulting from the leading-order coupling of the elec- 
tronic and nuclear motions, (ii) the relativistic corrections to 
the Schrodinger equation, (iii) the quantum electrodynamics 
(QED) corrections, and (iv) the residual retardation correc- 
tion. The largest contribution to the dimer energy ( — 1.615 
mK) is +0.226 mK repulsively from the correction (ii), and 
the total contribution is +0.103 mK with a mutual can- 
cellation among (i)-(iv) 11411 . The potential is referred to 
as "PCKLJS" (an acronym for Przybytek-Cencek-Komasa- 
Lach-Jeziorski-Szalewicz, the authors of Ref. IU4ID in a sub- 
sequent paper [21]; hereafter we use this acronym. 

The first purpose of the present work is to calculate, us- 
ing the potential PCKLJS, the binding energies of the trimer 
and tetramer ground and excited states, B^ ,B^ ,B 4 and 

B\ , respectively, together with the estimation of the individ- 
ual contributions of the corrections (i)-(iv). 

The large scattering length of 4 He- 4 He potential leads to 
universal properties in the four-body problem. An exam- 



pie is existence of the correlations between the different ob- 
servables. Thus, the second purpose of this work is to cal- 
culate the binding energies B^ ,B^ ,B\ and B A using 
various realistic 4 He potentials and investigate six types of 
the correlations between any two of the four energies. The 
potentials employed are PCKLJS and six other potentials: 
LM2M2 01TTY JH, HFD-B & HFD-B3-FCI1 fUJIl, 
S APT96 S^lli and CCSAPT07|S (see Ref. JH for a re- 
view of the recent study of the 4 He potential); in the last three, 
we choose the cases of the retardation corrections included. 

Recently, universal correlations between observables have 
been studied extensively in four-boson systems interacting 
with large scattering length J51 l3lM36ll . As for the specific 
4 He tetramers, the universal scaling functions for the cor- 
relations between the trimer and tetramer binding energies 
were obtained by the leading-order effective theory 04j |3J] and 
compared with the energies calculated using realistic 4 He po- 
tentials. However, due to the scarce calculation of the 4 He 
tetramer excited-state binding energy B^ ' using the realis- 
tic pair potential at that time, the correlations associated with 
B\ remained unexplored (see Fig. 4 of Ref. H). In the 
present work, we provide precise systematic results on the six 
types of correlations and demonstrate that the correlations are 
all linear over the range of binding energies relevant to 4 He 
atoms. We compare the result with that given by the leading- 
order effective theory 0|5|] for 4 He atoms. 

This paper is organized as follows: In Sec. II, we briefly 
present our calculational method GEM IlSl— TlOfl . Calculated 
results using the PCKLJS potential are presented in Sec. Ill 
together with those for the post-BO corrections (i)-(iv). In 
Sec. IV, using various 4 He potentials, we calculate the trimer 
and tetramer ground- and excited-state binding energies and 
discuss the correlations between them in comparison with the 
universal scaling functions obtained by the leading-order ef- 
fective theory for the 4 He atom. A summary is given in Sec. V. 



II. METHOD 

We employ the same ab initio variational method GEM as 
in the previous work [6] to solve the ground and excited states 
of the 4 He trimer and tetramer. Here, we review the method 
for the case of the tetramer. 

We take two types of Jacobi coordinate sets, K-type and 
ff-type (Fig. [TJ. For the K-type, xi = r 2 — ri, yi = r3 — 
|(ri + r 2 ) and zi = r 4 - |(ri + r 2 + r 3 ) and cyclically 
for {xj, y^, z^; i = 2, ..., 12} by the symmetrization between 
the four particles. For the iJ-type, X13 = r 2 — ri, yi3 = 
r 4 - r 3 , and z 13 = ±(r 3 + r 4 ) - ±(ri + r 2 ) and cyclically 
for{xi,yj,Zi; i = 14, ...,18}. 

The total four-body wave function ^4 is to be obtained by 
solving the Schodinger equation 



(H -£)# 4 = 



with the Hamiltonian 



H 



2/ix 



-V 



2(i ; 



■V 



2/i, 



-V 



4 

E 

1=»<7 



(2.1) 



V( rij ), (2.2) 




! 



'73 



O 

Yl3 



16 



(J3 



i=13 
H-type (/' =13-18) 



FIG. 1 : K-type and //-type Jacobi coordinates for the He tetramer. 
Symmetrization of the four particles generates the sets i — 1, ..., 12 
(isT-type) and i = 13, ..., 18 (H-type). 



where [i x — \vn, /i y = |m and \x z = |m on the K- 
type coordinates, and fi x = fj, y = \m and \i z = m on the 
i7-type ones, m being the mass of the 4 He atom. We take 
^ = 12.11928 KA 2 lUl. V(r y ) is the two-body 4 He- 4 He 
potential as a function of the pair separation ry = Vj — r^ . 

The wave function $4 is expanded in terms of the sym- 
metrized L 2 -integrable A'-type and H-type four-body basis 
functions: 

* 4 = v 4£ ) $iT n; * ) + E A i H J* ( «r H) > ( 2 - 3 ) 
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in which (x,;, y,;, Zj) is the i-th set of Jacobi coordinates. It is 
of importance that $Q^ m ' and $aj m ' are constructed on 
the full 18 sets of Jacobi coordinates; this makes the function 
space of the basis quite wide. 

The eigenenergies E and amplitudes A aK (A a J) are de- 
termined by the Rayleigh-Ritz variational principle: 

( $^ m « I H - E I * 4 } = 0, (2.6) 

( $^ m;H) I H - E I * 4 ) = 0, (2.7) 

where ax = !,•••, ax ax an< ^ a H = I, ■■-,a'^ ax . These 
equations result in the generalized matrix eigenvalue problem 
[Eqs.(3.8)-(3.10)ofpaperl]. 

We describe the basis function $L K (^a H ) m me form 

_[y ix (X 4 )^„(y 4 )] A ^,(^)] jM , 

(i = 1, ..., 12) (2.8) 

[y ix (S 4 )^(y 4 )] A 1^2 4 )] jM , 

(i = 13, ..., 18) (2.9) 



where olk specifies the set 

an = {cos or sin, u, n x l x , i 



iy Uy, IL Z 



l z ,A,JM}, (2.10) 



which is the same for the components i = 1, ..., 12; and 
similarly for a#, for all i = 13, ..., 18. J is the total an- 
gular momentum and M is its z component. In this paper, 
we consider the tetramer bound states with J = 0. There- 
fore, the totally symmetric four-body wave function requires 
(i) l x — even, l y + l z = even and A = l z for the if -type basis 
and (ii) l x = even, l y = even and A = l z = even for the 
H-type basis. 

In Eqs. (2.8) and (2.9), the radial functions are assumed as 

©■) =*'■«-<■'-•' *{zt'z§w 

i>n y i y (y) = y lv e~ {y/yny)2 , ( 2 - 12 ) 

(p n .i.{z) = z l - e-W*~-? (2.13) 

with geometric sequences of the Gaussian ranges: 
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(n y = l,...,n™*), (2.15) 

(n, = l,..,nr x )- (2.16) 



residual retardation correction V m t(r). Contributions of the 
individual corrections are discussed in Sec.IIIA. 

Use of PCKLJS for the dimer [ 14] gives the binding energy 
B 2 = 1.62 ± 0.03 mK, the average separation (r) = 47.1 ± 
0.5 A and the s-wave scattering length a = 90.42 ± 0.92 A. 
Experimental values of the quantities were reported 113911 as 
B 2 = l.l+oJ mK, (r) = 52 ± 4 A and a = 104±f 8 A, but 
the B 2 and a were calculated fl39ll from the observed value of 
(r) using rather crude models: B 2 — Ti 2 / (4m(r) 2 ) and a = 
2(r), where m is mass of 4 He atom. Much better estimates 
of what should be the values of B 2 and a correspon ding to 
the experimental (r) were recently obtained in Ref. [40], a 
follow-up paper to Ref. Q, to be B 2 = Ot^g mK and 
a = 100. 2ly g A, which are substantially close to and nearly 



consistent with the ab initio results 111411 mentioned above. 

Calculated binding energies of the trimer and tetramer for 
the PCKLJS potential are Bf ] = 131.84 mK, B { 3 1] = 2.6502 
mK (1.03 mK below the dimer), and Bf ] = 573.90 mK, 
B4 = 132.70 mK (0.86 mK below the trimer ground state). 
Some of the mean values of the trimer (tetramer) ground 
and excited states, as well as the binding energies mentioned 
above are summarized in Table HI (Table Hill . The PCKLJS po- 
tential gives slightly deeper binding of the trimer and tetramer 
than the LM2M2 potential does if. 



It shoud be emphasized that the GEM few-body calcula- 
tions need neither the introduction of any a priori pair corre- 
lation function (such as the Jastrow function) nor separation 
of the coordinate space into x < r c and x > r c , with r c be- 
ing the radius of a strongly repulsive core potential. Proper 
short-range correlation and asymptotic behavior of the to- 
tal wave function are automatically obtained by solving the 
Schrodinger equation (2.1) using the above basis functions for 
ab initio calculations. 

We take the same three- and four-body Gaussian basis func- 
tions as those employed in paper I. The numbers of the total 
bases are 4400 for the trimer and 29056 for the tetramer; the 
bases range from very compact to very diffuse with the Gaus- 
sian ranges in geometric sequences. 



III. THE PCKLJS POTENTIAL AND 4 He TRIMER AND 
TETRAMER 

The currently most accurate ab initio potential, the 
PCKLJS 1 14] potential, is given as a function of the 4 He pair 
separation distance r by 



A. Spatial structure of the tetramer 

We discuss the spatial structure of the tetramer excited state. 
In the study of four-boson states and their connection to the 
Efimov physics, the authors of Refs. IBU [3211 predicted that 
below each Efimov trimer a pair of tetramer states ( J 71 = + ) 
should exist and that the shallower member of the lowest-lying 
pair is dominantly composed of the ground-state trimer and a 
distant atom. It is shown, in the calculations by Lazauskas and 
Carbonell [42] and by the present authors [6] using the real- 
istic LM2M2 potential, that the above prediction is realized 
in the two bound states of the 4 He tetramer below the trimer 
ground state. 

The structure of the 4 He tetramer excited state is seen es- 
sentially in Fig. [2]for the overlap function 0± (z) between a 

tetramer state W^ (v = 0, 1) and the trimer ground state ^3 
which is defined as a function of the distance z between the 
trimer and the fourth atom: 



V(r) = V BO {r) + KdM + V Icl (r) + V QEB (r), (3.1) 

which are composed of the nonrelativistic BO potential ( Vbo) 
and the leading order coupling of the electronic and nuclear 
motions, that is, the adiabatic correction (V a d), relativistic 
corrections (V Te \), and quantum electrodynamics corrections 
(Vqed)- Besides them the Casimir-Polder retardation ef- 
fect [38], denoted as V ro t(r), can be optionally added to V(r). 
By the PCKLJS potential we mean the full V(r) plus the 
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Figure|2]indicates that the fourth atom is located in the trimer 
core region in the tetramer ground state but is far from the 
trimer in the excited state. This is also understood from the 
fact that, in Tables U and [TTJ the binding energy and the quan- 
tities (T) and (V) for the tetramer excited state are very close 
to those for the trimer ground state. 



TABLE I: The binding energies Bg'(v — 0, 1) and mean values 
of the 4 He trimer ground and excited states using the PCKLJS poten- 
tial 11411 including all the corrections. Tij stands for interparticle dis- 
tance and rio is the distance of a particle from the center of mass of 
the trimer. Cg is the asymptotic normalization coefficient defined 
by Eq.(2.25) in paper I. The conversion constant — = 12.11928 
KA 2 is taken. 



TABLE II: The binding energies B^ (v — 0, 1) and mean values of 
the 4 He tetramer ground and excited states using the PCKLJS poten- 
tial I14J1 including all the corrections. rij stands for interparticle dis- 
tance and TiG is the distance of a particle from the center-of-mass of 
the tetramer. CJ is the asymptotic normalization coefficient defined 
by Eq.(3.22) in paper I. The conversion constant — = 12.11928 
KA 2 is taken. 
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B^ (mK) 


131.84 
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FIG. 2: Overlap function O^ (z) in Eq. (3.2) between the tetramer 
state (v — 0, 1) and the trimer ground state as a function of the atom- 
trimer distance z. 



B. Relativistic and QED corrections 

The first four columns of Table Hill list the calculated dimer 
binding energy and the average interparticle distance at each 
level of theory [PCKLJS (a) to (h)], showing the contribu- 
tions of VbO) Kid) V le \, Vqed and the retardation corrections 
(denoted as "r.c"), which are different at different levels of 
theory [48]. The numbers in the first and second columns, 
given by Ref. O, are precisely reproduced by our calcula- 
tion. 

Using the potentials PCKLJS (a) to (h), we calculated the 
binding energies of the ground and excited states of the 4 He 



(0) r(1) r(0) 



,(D 



trimer and tetramer, B^ , -B 3 ,B\ and B\ , which are 
listed in Table III. In the tetramer ground-(excited-)state en- 



ergy, — B\ (—B\ ), the contribution from each correction is 
as follows: The retardation correction is +6.87 (+2.55) mK 
repulsively for VqOi but the residual (remaining) correction 
is only +0.16 (+0.07) mK for V BO + Kd + V+i + Vqed. 
Comparing (a) and (c), we see that the nonadiabatic correction 
(Kd) is -4.13 (-1.52) mK attractively. From (c) and (e), the 
relativistic correction (V Te \) is known as +9.37 (+3.48) mK. 
The QED correction (Vqed) is -1.20 (-0.46) mK from (e) 
and (g); the entire correction amounts to +4.20 (+1.57) mK. 

We remark that each correction for the tetramer excited- 
state energy (— B\ ') is approximately the same as the cor- 
responding correction for the trimer ground-state energy 
(— i?3 ). For example, the difference between e) and g), 
the QED correction, is -0.46 (-0.46) mK and that be- 
tween a) and h), the full correction, is +1.57 (+1.59) mK 
for — B\ (— B3 ). This is quite reasonable since the same 
explanation in the paragraph below Eq. (3.2) is applica- 
ble. A similar tendency is seen in the comparison between 
the correction for the trimer excited-state energy (— B a ) 
and that for the dimer (— B-z); for example, the QED cor- 
rection is —0.036 (—0.030) mK and the full correction is 

+0.12 (+0.10) mK for -Bf\-B 2 ). 



IV. UNIVERSALITY IN 4 He TRIMER AND TETRAMER 

Universal correlations between observables have been stud- 
ied systematically in four-boson systems interacting with 
large scattering length Il5l I31K36I1 . In this section we inves- 
tigate the correlations between the ground- and excited-state 
binding energies of the 4 He trimer and tetramer. We cal- 
culate the energies using various realistic 4 He- 4 He interac- 



TABLE III: Calculated binding energies of the He dimer, trimer and tetramer using the PCKLJS potential IU4I1 to demonstrate the contributions 
of Vbo , Vkd , Vrei , Vqed and the retardation correction, denoted as "r.c", appropriate for a given level of theory. PCKLJS-h is the full PCKLJS 
potential. B 2 and (r) are the binding energy and the average separation of the dimer, respectively. Calculated results for the dimer by Ref. IU4II 

are shown in the first and second columns. The conversion constant — = 12.11928 KA 2 is taken. 





Dimer [14] 
B 2 (r) 


Dimer 


Trimer 


Tetramer 


PCKLJS potential [14J 


B 2 


(r) 


Bf 


B?) 


bT 


B? 




(mK) 


(A) 


(mK) 


(A) 


(mK) 


(mK) 


(mK) 


(mK) 


(a) Vbo 


1.718 


45.77 


1.7181 


45.77 


133.43 


2.7724 


578.10 


134.27 


(b) Vbo + r.c. 


1.555 


47.92 


1.5549 


47.92 


130.85 


2.5776 


571.23 


131.72 


(c) Vbo + V ad 


1.816 


44.62 


1.8160 


44.62 


134.96 


2.8881 


582.23 


135.79 


(d) Vbo + V ad + r.c. 


1.648 


46.65 


1.6482 


46.65 


132.37 


2.6894 


575.33 


133.23 


(e) Vbo + V ad + V re i 


1.590 


47.43 


1.5896 


47.43 


131.44 


2.6194 


572.86 


132.31 


(f) Vbo + V ad + V rc i + r.c. 


1.611 


47.15 


1.6105 


47.15 


131.76 


2.6444 


573.69 


132.62 


(g) Vbo + V ad + Vol + Vqed 


1.620 


47.02 


1.6200 


47.02 


131.90 


2.6559 


574.06 


132.77 


(h) Vbo + V ad + V ro i + Vqed + r.c. 


1.615 


47.09 


1.6154 


47.09 


131.84 


2.6502 


573.90 


132.70 



TABLE IV: The binding energy B 2 and the average interparticle 
distance (r) of the dimer calculated using the seven He potentials. 
The conversion constant — = 12.11928 KA 2 is taken. The val- 

m 

ues reported in the literature are also shown, but the numbers in the 
parentheses were obtained by using — — 12.12 KA 2 and those in 
the square brackets were given with the use of the 4 He nuclear mass 
for m. The potential names are arranged in the increasing order of 
B 2 . 





This work 


Other work 


Potential 


B 2 


(r) 


B 2 


(r) 




(mK) 


(A) 


(mK) 


(A) Ref. 


LM2M2 


1.3094 


51.87 


(1.3035) 


(52.00) [6J 


TTY 


1.3156 


51.76 


(1.3096) 


(51.89) [411, 


HFD-B3-FCI1 


1.4475 


49.52 


1.448 


49.52 [2LJ 


CCSAPT07 


1.5643 


47.78 


1.56 


47.8 [29] 


PCKLJS 


1.6154 


47.09 


1.615 


47.09 [14] 


HFD-B 


1.6921 


46.07 


(1.6854) 


(46.18) [4LJ 


SAPT96 


1.7443 


45.45 


[1.713] 


[45.8] [28J 



tions which include the PCKLJS potential and six other po- 
tentials LM2M2, TTY, HFD-B, HFD-B3-FCI1, SAPT96 and 
CCSAPT07 mentioned in Sec. I; in the last three, we choose 
the cases in which the retardation corrections are included. 

We first calculated, using the seven potentials, the binding 
energy B2 and the average interparticle distance (r) of the 
dimer and listed them in Table [IV] together with the values 
reported in the literature. The scattering lengths are not listed, 
but they range between 87.92 A [26] for SAPT96 and 100.23 



TABLE V: The binding energies of 4 He trimer and tetramer ground 
and excited states calculated with the use of the seven 4 He potentials. 
The conversion constant — = 12.11928 KA 2 is taken. 



Trimer 



Tetramer 



Potential 



B 



(0) 



B 



W 



B 



(0) 



B 



CD 



(mK) (mK) 



(mK) (mK) 



LM2M2 


126.50 


2.2779 


559.22 


127.42 


TTY 


126.45 


2.2844 


558.70 


127.37 


HFD-B3-FCI1 


129.00 


2.4475 


566.12 


129.89 


CCSAPT07 


131.01 


2.5890 


571.67 


131.88 


PCKLJS 


131.84 


2.6502 


573.90 


132.70 


HFD-B 


133.08 


2.7420 


577.34 


133.94 


SAPT96 


134.02 


2.8045 


580.01 


134.86 



A lU for LM2M2 (90.42 A UM for PCKLJS). Here, the 
names of the seven potentials are arranged from the top to the 
bottom in the increasing order of i?2- 

The 4 He trimer and tetramer ground- and excited-state 
binding energies, B 3 ,B 3 ,B^ and B\ , are calculated 
with those potentials and are listed in Table V. The val- 
ues of each binding energy appear in the increasing order 



(o) p (o) 



?(i) 



as B2 does in Table IV except for B 3 ,B\ and B\ for 
LM2M2 and TTY. This exception is reasonable because TTY 
is slightly more attractive for r > 2.65 A than LM2M2, but 
slightly more repulsive for r < 2.65 A; namely, it is possi- 
ble that TTY generates larger binding energies than LM2M2 
in loosely bound systems (the dimer and the trimer excited 
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FIG. 3: (Color online). The correlations between the ground- and excited-state binding energies of the He trimer and tetramer. (a) B\ -B, 



(1) 7,00 



(b)Bf 



-B\ , (c) Bg -B\ and (d) B^ -B\ correlations. The energies are normalized by the dimer energy B^- All the 14 data points 



are obtained by the present calculation for various 4 He potentials. The seven closed circles, from the right to the left, denote the results for 
the seven potentials in Table V from LM2M2 down to SAPT96, respectively; the one designated by an arrow is for PCKLJS. The seven open 
squares, from the left to the right, show the calculation for each level of the PCKLJS potential in Table I in the order of PCKLJS (c), (a), (d), 
(g), (f), (e), and (b), respectively. The dotted (red) linear line is the linear least squares fit to the 14 data points; see Eqs. (4. l)-(4.4). The 
dashed (blue) line in panel (b) is the prediction by the dimerlike-pair model, Eq. (4.8). The solid line, taken from Fig. 4 and Eqs. (39)-(42) in 
Ref. |5], shows the universal scaling curve obtained by the leading-order effective theory for the 4 He trimer and tetramer. 



states) but brings about smaller binding energies in compactly 
bound systems (the trimer and tetramer ground states and 
the tetramer excited state that is dominantly composed of the 
compact trimer ground state and a distant 4 He atom). We note 
that if all the binding energies are normalize by B^, they ap- 
pear in the increasing order in Table V. 



A. Linear correlations 

Correlations between the binding energies in three- and 
four-body systems were first observed in nuclear physics and 
are known as the Tjon line 114311 . which refers to the approxi- 
mately linear correlation between the binding energies of the 
triton and the a particle for various nucleon-nucleon poten- 
tials. Recently, the nuclear Tjon line was discussed in the 
context of the effective field theory of short-range interactions 



and low-momentum nucleon-nucleon potentials 1144 14511 . The 
Tjon lines for the 4 He trimer and tetramers were investigated 
in Refs. [J, |5|] over the range of binding energies relevant to 
4 He atoms on the basis of the leading-order effective theory. 
However, due to the scarce calculation of the 4 He tetramer 
excited-state binding energy B\ with the realistic 4 He po- 
tential at that time, the correlations associated with B\ ' re- 
mained unexplored. 

We consider all six kinds of the correlations between two of 



and B A , that are 



the four binding energies, B^ , -B 3 , B\ 
calculated using the seven 4 He potentials in Table V and the 
seven potentials of PCKLJS (a) to (g) in Table I. The binding 
energies are normalized by B2 , which is different for differ- 
ent 4 He potentials; this is due to the fact that the experimental 
value of the dimer binding energy has not been precisely ob- 
tained as mentioned at the beginning of Sec. III. 
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FIG. 4: (Color online), (a) B\ -B\ and (b) B\ -B\ correlations. Meanings of the data points (the present calculation) are the same as 
in Fig. [3] The dotted (red) linear line is the linear least squares fit to the 14 data points; see Eqs. (4.5) and (4.6). The solid line shows the 
universal scaling curve obtained by the leading-order effective theory for the 4 He atoms; we derived the line from Eqs. (39)-(42) of Ref. U 
(see footnote [49]). The dashed line in (a) is another universal scaling curve obtained in Ref. [4]. 



Figures Ha) to Hd) illustrate the B^-B^, B^-B^, 
i?3 -B± and B 3 -B 4 ' correlations, respectively. The 
14 data points are given by the present calculation for vari- 
ous 14 potentials mentioned above. The dotted (red) lines in 
Figs. |3 a) to[3jd) are obtained by the linear least squares fit- 
ting to the data points and are represented by the following 
equations, respectively: 



B 



(i) 



Bo 



151.4- 



B 



(i) 



Bo 



166.1, 70 < 



B 



(i) 



Bo 



<100, (4.1) 



B 



(i) 



Bo 



1.011- 



B 



(0) 



Bo 



B 



(0) 



Bo 



715.9- 



B 



(i) 



Bo 



B 



(0) 



Bo 



= 4.778- 



B 



(o) 



Bo 



0.3694, 



819.0, 



34.64, 



1.5 < -2- 
Bo 

^S« 
70 < 3 



Bo 



< 1.8, (4.2) 



< 100, (4.3) 



1-5 < -4- 

Bo 



< 1.8. (4.4) 



FiguresHa) and 4(b) plot the B^'-B^ 1 and Bf'-BY' cor 



iW-BP and bP-BP' 
relations. With the least squares method, the data points are 
fitted by the dotted (red) lines that are represented by 



B 



(i) 



Bo 



0.006679- 



B 



(0) 



B 



(i) 



Bo 



Ba 

0.2116 ' 



B 2 

(0) 



1.095, 



Bo 



6.961, 



R (0) 

300 <^±- 
Bo 



< 100, (4.5) 



<440. (4.6) 



In Figs. [3] and |4] the scattering of the data points about 
the fitted linear line is very small: Representing the data 
points by {(Xi, j/j), i — 1, ..., 14} and the fitted linear func- 
tion by y = f{x), we define relative deviation at each Xi by 
\yi — f(xi)\/yi- The average values of the relative deviation 
in Figs. 3(a), 3(b), ...., 4(b) are respectively 0.093%, 0.0032%, 



0.1 1%, 0.019%, 0.030% and 0.015%. We remark that, among 
Eqs. (4. 1 )-(4.6), any three equations can be reproduced by the 
other three (linearly dependent) with very small errors. This 
comes from the fact that the six kinds of correlations are all 
linear for various potentials. 

It is unexpected that all the calculated results (the data 
points) fall so strictly on a straight line over the range of bind- 
ing energies relevant for 4 He atoms; we emphasize that the re- 
sults are obtained by using different potentials, not by chang- 
ing parameter(s) in a specific potential. 

It is of interest to note that the slope of the dotted (red) line 
in Fig. 3(d) for the correlation between trimer and tetramer 
ground-state binding energies is 4.778 [seeEq. (4.4)], which is 
close to the slope of the nuclear Tjon line (« 5.0 114411 ) for the 
correlation between three- and four-nucleon binding energies 
using various nucleon-nucleon potentials. 

Another similarity between the 4 He tetramer and the 4 He 
nucleus is seen in the comparison of the overlap function 
O4 (z) in Fig. 2 ( 4 He tetramer) with that in Fig. 5 ( 4 He nu- 
cleus). The behavior of O4 (z) is quite resemble each other 
although the sizes of the systems are very different. The first 
excited + state of the 4 He nucleus is known to be composed 
of the three-nucleon core and a loosely coupled nucleon [46]. 
The observed cross sections of the electron inelastic scatter- 
ing, which drastically excites the compact ground state to the 
diffuse excited state, is well explained by the GEM four-body 
calculation by the present authors [46]. 

In the study of weakly bound four-boson states (not 
specifically for 4 He atoms) at the unitary limit, von Stecher et 



cd. 132J] obtained B^/B^' « 4.58 and B, 



? (D 



while Deltuva 

B^/B^ 



B 



(o) /n (o) 



gave tSl''/B 3 
1.0023, and Hadizadeh et al. 
4.6. In the 4 He atoms, 



/B^> a 1.01, 

» 4.611 and 

113311 reported 

we can estimate 



the ratio, approximately from Eqs. (4.4) and (4.2), 



as Bf'/Bf> = 4.778 - M.UB 2 /B$ 
B^/B^ = 1.011 - 0.3694 B 2 /B {0] 



(o) 



« 4.3^.4 and 
1.006-1.007 over 



the range of binding energies relevant for the 4 He atoms. 
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FIG. 6: Schematic picture of the dimerlike-pair model for the trimer 
and tetramer excited states in the asymptotic region. 



FIG. 5: Overlap function 0±{z) between the 4 He nucleus (v = 
0,1) and the three-nucleon ground state as a function of the distance 
z between the three-nucleon core and the fourth nucleon. Taken from 
Ref. I4ql . Note that this figure is quite similar to Fig. 2 for the 4 He 
tetramer. 



B. Universal scaling functions 

The solid lines in Fig. [3] illustrate the universal scaling 
functions relating the tetramer energies to the trimer ener- 
gies, which were calculated by the leading-order effective the- 
ory for the 4 He atoms; the lines are taken from Eqs. (39)- 
(42) and Fig. 4 in Ref. ||5[]. To obtain the energies, Platter 
et al. yt] constructed an effective 4 He- 4 He potential includ- 
ing both two- and three-body contact interactions. The two 
parameters of the effective potential were determined to re- 
produce the binding energy of the dimer ground state and the 
trimer excited state. They solved the three- and four-body 
Faddeev-Yakubovsky equations with the effective potential. 



>(0) n (o) 



(0) p(l) 



Although the B^ -B% and B\ -B\ ' correlations are not 
explicitly given in Ref. J5[], we derived the functions for the 
correlations IH] using Eqs. (39)-(42) in Ref. [5] and plot 
them in Fig.|4]with the solid lines. The dashed line in Fig.Ua) 
is another universal scaling curve for the 4 He trimer given in 
Fig. 2 of Ref. H. 

The solid lines in Figs. [2 a), Ob) and|4ja), associated with 



>(o) R (i) 



i(i) 



B\ , B3 ' and B\ ' , are close to the calculated data points. 
Therefore, origin of the non-negligible deviation of the solid 
lines in Figs. |3jc) and |3jd) and Fig.Hfb) from the data points 
is attributed to B\ ' given by the leading-order effective the- 
ory H. We estimate that the theory underestimates B\ by 
some 70-80 mK (~ 12-14% of B { ° ] ) while it overestimates 



jW 



(i), 



B\ ' by about 2-3 mK (~ 2-3% of B\ ') compared with the 
calculation using the realistic 4 He potentials. 



C. Dimerlike-pair model 

In Fig.[3b) for the B% -B A correlation, the dashed (blue) 
line, predicted by the dimerlike-pair model y6|], is close to the 
14 data points with almost the same quality as the dotted (red) 
line of the least squares fit. 

We briefly recapitulate the model. Firstly, for the trimer 
excited state in Fig. 6(a), the model indicates that (i) parti- 
cle a, located far from b and c (dimer), which are loosely 
bound, is little affected by the interaction between b and c; 
(ii) therefore, the pair a and b at a distance x is asymptotically 
dimerlike; (iii) since x ~ y asymptotically, the wave function 
of particle a along y, exp(— fcg y)/y, is the same as that of 
the dimer, exp(— k2x)/x; hence we have a relation fcg = k 2 - 
The binding wave numbers are related to the binding energies 

as 

k 2 = y/2(i x B 2 /h, 



fcW _ 

n 3 — 



,(i) 



B 2 )/n, 



2/j, v (B 3 

where /i x = \m and [i y — |to are the reduced mass asso- 
ciated with the coordinates x and y, respectively. Using the 
relation fc 3 = k 2 , the model predicts 



B, 



(i) 



Bo 



B 2 



-AS 



(i) 



B 2 



(4.7) 



where AS 3 = B^ ' — B 2 is the binding energy measured 
from the dimer. 

Similarly, the model predicts the tetramer excited-state en- 
ergy as follows: Asymptotically, in Fig. 6b, particle a decays 
from the trimer (6 + c + d) as exp(- 



-k\ z)/z with 



n/4 — 



2^ z {B { l ) ~B < i ) )/n, 



where fx z — |m is the reduced atom-trimer mass. Using 
ki = k 2 by the same reason as above, the model predicts 



B 



(i) 



B 



(0) 



Bo 



B 2 



/AB 



(i) 



V B 2 



(4.8) 



TABLE VI: Comparison of the trimer (tetramer) excited-state bind- 
ing energy B^ (B A ) between the dimerlike-pair-model prediction 

in the column "model" and the present three-body (four-body) cal- 

As (i) . „(« 

culation. The model predicts B 3 — 



| and 



Generally, for the 4 Hejv, the model suggests a relation 



B 



(i) 

JV 



D 



(0) 
AT-l 



Bo 



Bo 



N 



2(N-iy 



(4.9) 



Model 



Present calculation 



Potential 



B 



CD 



B 



(i) 



o(l) 
n 3 



o(l) 



(mK) (mK) (mK) (mK) 



LM2M2 2.29 127.37 2.2779 127.42 0.74 0.70 

TTY 2.30 127.33 2.2844 127.37 0.74 0.70 

HFD-B3-FCI1 2.53 129.96 2.4475 129.89 0.69 0.61 

CCSAPT07 2.74 132.05 2.5890 131.88 0.66 0.56 

PCKLJS 2.83 132.91 2.6502 132.70 0.64 0.53 

HFD-B 2.96 134.21 2.7420 133.94 0.62 0.51 

SAPT96 3.05 135.18 2.8045 134.86 0.61 0.48 



where AB\ = B\ — B$ is the binding energy measured 
from the trimer ground state. 

We emphasize that the relations (4.7) and (4.8) are interac- 
tion independent. Our assumption is only that the interparti- 
cle distance in the dimer (trimer) is larger than the interaction 
range, which is fulfilled in the present case. The relation (4.8) 
is plotted in Fig.[3jb) by the dashed (blue) line, which almost 
overlaps with the dotted (red) line for the linear fit to the 14 
data points by the present few-body calculation. 

The trimer and tetramer excited-state binding energies 
predicted by the dimerlike-pair model are summerized in 
Table VI in comparison with the results of the present three- 
and four-body calculations using various 4 He interaction. Er- 
ror of the model prediction is known to be 0.01-0.25 mK in 
B^iAB^) and 0.05-0.32 mK in b[ 1 \Ab[ 1) ) compared 
with the three-(four-)body calculation. The model works well 
when we compare, in Fig. 3(b) (the B 3 -B\ correlation), 
the 14 data points with the dashed (blue) line and with the 
solid line. 

However, strictly speaking about the model, it predicts 
neither B ( 3 \b^ ] ) nor the threshold energy B 2 (Bf ] ) but 



where the last term is ratio of the reduced mass of the dimer- 
like pair ( |m) to that of the 4 He- 4 Hejv-i system {^^-m). 



V. SUMMARY 

Using the Gaussian expansion method for ab initio varia- 
tional calculations of few-body systems [8-10], we have cal- 
culated the binding energies of the 4 He trimer and tetramer 



(o) R (i) R (o) 



?d) 



ground and excited states, B 3 ,B 3 ,B\ ' and B\ ' , with the 
use of the currently most accurate 4 He potential proposed by 
Przybytek et al. [14], called the PCKLJS potential. This is an 
extension of our previous work [6] using the LM2M2 poten- 
tial. Employing the PCKLJS, LM2M2, TTY, HFD-B, HFD- 
B3-FCI1, SAPT96 and CCSAPT07 potentials, we have cal- 
culated the four kinds of the binding energies and investigated 
the correlations between any two of them. 

The main conclusions are summarized as follows: 



(i) We obtained, using PCKLJS, B 



(0) 



131.84 mK, 



B 



(i) _ 



2.6502 mK (1.03 mK below the dimer), B 



(o) _ 



573.90 mK and B\ L ' = 132.70 mK (0.86 mK below the trimer 
ground state). This potential includes the adiabatic, relativis- 
tic, QED and residual retardation corrections. Contributions 



(0), 



>dh 



of the corrections to the tetramer eigenenergy — B\ {—B\ ') 
are, respectively, -4.13 (-1.52) mK, +9.37 (+3.48) mK, 
-1.20 (-0.46) mK and +0.16 (+0.07) mK; the entire cor- 
rection is +4.20 (+1.57) mK. 

(ii) The six correlations B { 3 ] -B^ ] , B^-B^, B { ^-Bf\ 

Bf ] -Bf\ B^-B^ and B^-B^K are observed to be all 
linear for the various 4 He potentials mentioned above over the 
range of binding energies relevant for the 4 He atoms. They 
may be called the generalized Tjon lines. The universal scal- 
ing curves given by the leading-order effective theory for the 
4 He atoms [5] locate closely to the presently-obtained linear 



»(i)/ 



>(i) 



lines in the B^-B^ \ Bk ' B\ ' ' and U 



i (1) - 



i 1] 



(0) o(l) 



B, ; correlations, 

(o) D (o) 



does predict their difference AB y 3 '(AB\ '). Therefore, the but deviate non-negligibly from the lines in the B^'-B, 
error or accuracy of the model prediction might be judged 



»(!)/ 



,(D 



M -B& an d B\ -B A correlations. The latter deviations 



also on the basis of AB y 3 (AB\ ) . For this purpose, one can are attributed to the fact that the leading-order effective theory 

A R^ 1 ) A R^ 1 ) (Q) (0) 

compare in Table VI the value of =ga_ (^-) by the three- underestimates B\ ' by some 70-80 mK (-12-14% of B\ ') 



(four-)body calculation with the value, | (|), by the model. 



ABi 



-AS 



The relative error of the model prediction for B 3 ( B 4 ) 
is estimated as 1-19% (5-28%)|H3]. This is not very small, 
but we remark that the simple model provides, without any 

AB (1) 

elaborate calculation, a reason why the quantity B 3 ap- 



compared with the present calculations using the realistic 4 He 
potentials. 

(iii) As long as the binding energies of the excited states 
of the trimer and tetramer, the interaction-independent predic- 

tion, B 



I md — 2 



pears close to | (0.61-0.74 in Table VI), and 
to | (0.48-0.70). 



AB 



is close 



4, by the dimerlike-pair 

model [6] explains the B2-S3 andi?^ -B\ correlations for 
the various 4 He potentials, with B 2 being the dimer binding 
energy. 



10 



Acknowledgement 

We thank Krzysztof Szalewicz for providing us with 
codes producing the 4 He potential (PCKLJS) used in this 



work. The numerical calculations were performed on 
HITACHI SRI 6000 at KEK and YIFP. 



[1] V. Efimov, Yad. Fiz. 12, 1080 (1970) [Sov. J. Nucl. Phys. 12, 

589 (1971)]; Nucl. Phys. A 210, 157 (1973). 
[2] V. Efimov, Few-Body Syst. 51, 79 (201 1). 
[3] E. Braaten and H-W. Hammer, Phys. Reports, 428, 259 (2006). 
[4] E. Braaten and H-W. Hammer, Phys. Rev. A 67, 042706 (2003). 
[5] L. Platter, H.-W. Hammer and Ulf-G. Meissner, Phys. Rev. A 

70, 052101 (2004). 
[6] E. Hiyama and M. Kamimura, Phys. Rev. A 85, 022502 (2012). 
[7] R.A. Aziz and M.J. Slaman, J. Chem. Phys. 94, 8047 (1991). 
[8] M. Kamimura, Phys. Rev. A38, 621 (1988). 
[9] H. Kameyama, M. Kamimura and Y. Fukushima, Phys. Rev. C 

40,974(1989). 
[10] E. Hiyama, Y Kino and M. Kamimura, Prog. Part. Nucl. Phys. 

51, 223 (2003). 
[11] E. Hiyama and T. Yamada, Prog. Part. Nucl. Phys. 63, 339 

(2009). 
[12] E. Hiyama et ah, Prog. Theor. Phys. Supplement 185, 106 

(2010); 185, 152 (2010). 
[13] E. Hiyama, Few-Body Systems 53 (2012), Online First, 14 

March, 2012. 
[14] M. Przybytek, W. Cencek, J. Komasa, G. Lach, B. 

Jeziorski and K. Szalewicz, Phys. Rev. Lett. 104, 

183003 (2010); see also the supplementary material at 

http://link.aps.org/supplemental/10.1103/PhysRevLett.104.183003 

for extended tables and the parameters of the fit. 
[15] J. Fischer and B. Fellmuth, Rep. Prog. Phys. 68, 1043 (2005). 
[16] L. Pitre, M. R. Moldover, and W. L. Tew, Metrologia 43, 142 

(2006). 
[17] J. B. Mehl, Comptes Rendus Physique 10, 859 (2009). 
[18] S. Goyal, D. L. Schutt, and G. Scoles, Phys. Rev. Lett. 69, 933 

(1992). 
[19] M. Hartmann, R. E. Miller, J. P. Toennies, and A. Vilesov, Phys. 

Rev. Lett. 75, 1566 (1995). 
[20] S. Y Chang and M. Boninsegni, J. Chem. Phys. 115, 2629 

(2001). 
[21] V. Spirko, S.P.A. Sauer and K. Szalewicz, "On the relation be- 
tween properties of helium dimer bound state", preprint (201 1). 
[22] K.T. Tang, J.P Toennies and C.L. Yiu, Phys. Rev. Lett. 74, 1546 

(1995). 
[23] R.A. Aziz, F.R.W. McCourt and C.C.K. Wong, Mol. Chem. 

Phys. 61, 1487 (1987). 
[24] T van Mourik and J. H. van Lenthe, J. Chem. Phys. 102, 7479 

(1995). 
[25] R. A. Aziz, A. R. Janzen, and M. R. Moldover, Phys. Rev. Lett. 

74, 1586 (1995). 
[26] A.R. Jansen and R.A. Aziz, J. Chem. Phys. 107, 914 (1997). 
[27] H. L. Williams, T. Korona, R. Bukowski, B. Jeziorski, and K. 

Szalewicz, Chem. Phys. Lett. 262, 431 (1996). 
[28] T. Korona, H. L. Williams, R. Bukowski, B. Jeziorski, and K. 

Szalewicz, J. Chem. Phys. 106, 5109 (1997). 
[29] M. Jeziorska, W. Cencek, K. Patkowski, B. Jeziorski, and K. 



Szalewicz, J. Chem. Phys. 127, 124303 (2007). 
[30] K. Szalewicz, Int. Rev. Phys. Chem. 27, 273 (2008). 
[31] H.-W. Hammer and L. Platter, Eur. Phys. J. A 32, 1 13 (2007). 
[32] J. von Stecher, J.P. DTncao and C.H. Greene, Nature Phys., 5, 

417 (2009). 
[33] J.P. DTncao, J. von Stecher and C.H. Greene, Phys. Rev. Lett. 

103, 033004 (2009). 
[34] M.T. Yamashita, L. Tomio, A. Delfino and T. Frederico, Euro- 

phys. Lett., 75, 555 (2006). 
[35] M.R. Hadizadeh, M.T. Yamashita, L. Tomio, A. Delfino and T 

Frederico, Phys. Rev. Lett., 107, 135304 (2011). 
[36] T Frederico, L. Tomio, A. Delfino, M.R. Hadizadeh and M.T. 

Yamashita, Few-Body Syst., 51, 87 (2011). 
[37] V. Roudnev and M. Cavagnero, J. Phys. B: At. Mol. Opt. Phys. 

45,025101(2012). 
[38] H.B.G. Casimir and D. Polder, Phys. Rev. 73, 360 (1948). 
[39] R. Grisenti etal, Phys. Rev. Lett. 85, 2284 (2000). 
[40] W. Cencek, M. Przybytek, J. Komasa, J.B. Mehl, B. Jeziorski 

and K. Szalewicz, "Effects of adiabatic, relativistic, and quan- 
tum electrodynamics interactions in helium dimer on ther- 

mophysical properties of helium", submitted to J. Chem. 

Phys.(2012) (a follow up paper of Ref. 03). 
[41] E.A. Kolganova, A.K. Motovilov and W. Sandhas, Few-Body 

Syst. 51, 249(2011). 
[42] R. Lazauskas and J. Carbonell, Phys. Rev. A 73, 062717 (2006). 
[43] J.A. Tjon, Phys. Lett. B 56, 217 (1975). 
[44] A. Nogga, S.K. Bogner and A. Schwenk, Phys. Rev. C 70, 

061002 (2004). 
[45] L. Platter, H.-W. Hammer and U.-G. Meissner, Phys. Lett. B 

607, 254 (2005). 
[46] E. Hiyama, B.F Gibson and M. Kamimura, Phys. Rev. C 70, 

031001(R) (2004). 
[47] A. Deltuva, Phys. Rev. A 82, 040701(R) (2010). 
[48] For example, the retardation correction for Vbo is 0. 16 mK, but 

that for Vbo + Vad + Koi + Vqed is 0.005 mK in the dimer; 

namely, the residual contribution becomes much smaller in the 

latter. See Ref. Q3] for the details. 
[49] In Ref. [J, one can obtain an equation to relate B\ and 

B 3 by eliminating B 4 from Eqs. (39) and (40) and one 

more equation from Eqs. (41) and (42) by eliminating B A . 

Since the resultant two equations are slightly different from 

each other in the coefficients, we averaged them and obtained 



b 2 



B 2 



0.2504 -|— + 5.480. They are plotted in Fig.|4]with solid lines 
[50] One sees in Table VI that the error of the model increases with 
increasing B2 from the top to the bottom. However, investiga- 
tion of the reason for this behavior is out of the scope of the 
present work. 



